Black hole initial data on hyperboloidal slices 



o 
o 

(N 
> 

o 
O 
1^ 



cr: 



> 

m 

cn 
o 

OS 
O 



X 



Luisa T. Buchman/'^ Harald P. PfeifFer,^'^ and James M. Bardeen^ 

^Theoretical Astrophysics, California Institute of Technology, Pasadena, CA 91125 
^Center for Relativity, University of Texas at Austin, Austin, Texas 78712 
'^Canadian Institute for Theoretical Astrophysics, University of Toronto, Toronto, Ontario MBS 3H8, Canada 
'''Physics Department, University of Washington, Seattle, WA 98195 
(Dated: October 7, 2009) 

We generalize Bowen-York black hole initial data to hyperboloidal constant mean curvature slices 
which extend to future null infinity. We solve this initial value problem numerically for several cases, 
including unequal mass binary black holes with spins and boosts. The singularity at null infinity 
in the Hamiltonian constraint associated with a constant mean curvature hypersurface does not 
pose any particular difficulties. The inner boundaries of our slices are minimal surfaces. Trumpet 
configurations are explored both analytically and numerically. 
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I. INTRODUCTION 

For asymptotically flat spacctimcs, gravitational radi- 
ation is well-defined only at future null infinity {.f^) 
P, 0, 0|- Consequently, numerical relativists extract- 
ing gravitational waves from binary black hole simula- 
tions ideally would like to include in their compu- 
tational domains, so that the Bondi news function [l|, |j| 
(which contains the gravitational wave information) can 
be computed. Additionally, extending the simulation 
to would make it unnecessary to deal with grav- 
itational wave extraction at a finite distance, or with 
artificial outer boundaries on a truncated domain, two 
very complicated aspects of black hole simulations (see, 

e.g. s i a i, i, 113, El, m, El, 13 ) . 

Null infinity can be included in the computational do- 
main via a compactificd radial coordinate on the null 
hypersurfaces of the characteristic initial value prob- 
lem [l^ , but the null hypersurfaces are subject to caus- 
tic singularities, particularly in the strong fields around 
black holes. The caustics can be avoided by Cauchy- 
characteristic matching, but an appealing alternative is 
solving the Cauchy problem on conformally compactified 
hyperboloidal spacelike slices. These behave like conven- 
tional 3-1-1 slicing in the vicinity of the sources, but 
smoothly become asymptotically null as they approach 
null infinity at a finite coordinate distance |16l . [ItI . [ist . 
Fricdrich jl7| derived a system of symmetric hyperbolic 
evolution equations based on the Bianchi identities for 
the conformal Weyl tensor which are regular at null in- 
finity on hyperboloidal hypersurfaces provided certain 
smoothness conditions are satisfied. This system has 
been used with some, but limited, success in numeri- 
cal calculations (see, e.g. [3). More recently, evolu- 
tion schemes have been proposed which directly evolve 
the conformal metric and the conformal factor through 
the conformally compactified the Einstein equations on 
hyperboloidal hypersurfaces 2^21|. They are not man- 
ifestly regular at , but in [2ll|. Moncrief and Rinne 
derive regularity conditions based on the constraint equa- 
tions which deal successfully with the singularity of the 



conformal factor at and can be imposed in a numer- 
ical implementation. 

Numerical evolution schemes on hyperboloidal hyper- 
surfaces extending to J^'^ require initial data. The ini- 
tial value problem on asymptotically null hyperboloidal 
constant mean curvature (CMC) slices turns out to be 
remarkably similar to the corresponding problem on 
asymptotically fiat zero mean curvature slices, and can 
be attacked similarly to the well-known conformally flat 
Bowen-York [2^ initial data. (The Bowen-York initial 
data presented in [23] also forms the basis of puncture 
data [23]^). Therefore, in this paper, we shall refer to 
the initial data constructed as hyperboloidal Bowen- York 
data. 

There are four main points to this paper: (i) to 
lay down the formalism for constructing hyperboloidal 
Bowen-York black hole initial data on CMC slices con- 
taining one or more black holes with arbitrary masses, 
spins and boosts, (ii) to give rules for choosing the var- 
ious free parameters entering the formalism, (iii) to un- 
derstand the physical meaning of the free parameters 
and (iv) to discuss the physical interpretation of the con- 
structed solutions, which is different for CMC slices than 
for traditional maximal slices. 

Part of the initial data construction is the numerical 
solution of an elliptic equation arising from the Hamilto- 
nian constraint. On hyperboloidal slices, this equation is 
formally singular at the outer boundary Neverthe- 
less, the employed spectral elliptic solver [2J| does not 
exhibit any problems while computing the solution. 

Our initial data is formulated with a singularity- 
avoiding minimal surface boundary condition, on an 
Einstein- Rosen bridge connecting two asymptotically flat 



^ Puncture data differs from inversion symmetric data in the han- 
dling of the black hole singularities; we do not discuss the punc- 
ture treatment of singularities on CMC slices, although all our 
results concerning the constraint equations will carry over to such 
a treatment. 
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ends. In the limit in which the conformal radius ap- 
proaches zero, a "trumpet" configuration is formed in 
which the Einstein-Rosen bridge becomes infinitely long 
in proper distance. Hyperboloidal slicings that contain 
trumpets have been examined in Refs. [251. [26t. Trumpet 
initial data is of interest to numerical relativists using 
the moving puncture approach [2^ . [28j because puncture 
initial data evolve quickly toward a trumpet configura- 
tion [2i,[33,[3l|. 

The organization of this paper is as follows: Sec- 
tion|TT]presents and analyzes the formalism for construct- 
ing hyperboloidal Bowen-York initial data. Specifically, 
Sec. Ill Al presents the initial value formalism on CMC 
slices and Sec. Ill Bl gives the particulars in the special 
case of a Schwarzschild black hole. Sec. Ill CI discusses 
the parameter space yielding a minimal surface at the 
inner boundary of a Schwarzschild black hole in CMC 
slicing. Sec. Ill Dl presents the construction of initial data 
for single black holes, with and without Bowen-York spin 
and boost parameters. The physical interpretation of the 
Bowen-York parameters is also discussed. In Sec. Ill El 
the methods presented for single black holes are general- 
ized to multiple black holes. Sec. Ill Fl details the solution 
to the Hamiltonian constraint equation when the inner 
boundary is a trumpet. Sec. IIIII presents our numeri- 
cal solutions for hyperboloidal Bowen-York initial data. 
Specifically, Sees. IIII Al IIIIBI and IIII CI give solutions for 
single black holes that are spherically symmetric, spin- 
ning, or boosted, respectively. Sec. IIIIDI presents nu- 
merical results for two unequal mass black holes with 
arbitrarily oriented spins and boosts. Finally, Sec. IIVI 
discusses our results. An Appendix is included, in which 
conditions are derived for inversion symmetry on both 
CMC and maximal slices. 



II. ANALYTICS 
A. Initial value formalism on CMC slices 

With a standard 3+1 decomposition [s^, H^, the 
spacetime metric is written as 

ds^ = -a^ dt^ + gij {dx' + 13' dt) {dx^ + dt), (1) 

where gij represents the induced metric on the t = const 
hypersurface Et, a the lapse- function and /3* the shift- 
vector. The extrinsic curvature is defined as 

Ktra = 2'^'<')ri 9hn, (2) 

where £(4)„ is the Lie derivative along the hypersurface 
unit-normal '^'n^. In Eq. ([2|), as throughout this paper, 
we employ the sign convention of Wald |34l | , resulting in a 
positive mean curvature K = gijK'^^ for the cases consid- 
ered. This sign-convention differs from the conventions 
of Misner, Thorne & Wheeler [s^, which are more com- 
monly used in numerical relativity (as. for instance, in 
Ref. li]). 



Einstein's vacuum constraint equations are 

n + K^- K,j K'^ = 0, (3) 

and 

W,{K'^ ^g'^K)^0, (4) 

where Vi is the covariant derivative with respect to g^j, 
and TZ is the Ricci scalar associated with . 

We now perform a conformal transformation which 
plays a dual role. On one hand, it allows a conformal 
compactification (Penrose [IJ), placing ^+ at a finite 
value of a compactificd radial coordinate, absorbing the 
resulting metric singularities into the conformal factor. 
On the other hand, it allows the Einstein constraints to 
be recast as elliptic equations following the standard pro- 
cedure of the conformal method [ssl . l3a | . 

The conformal metric gij and conformal factor are 
given by 

g,,^^~^~g,„ g'^=VL^~g'K (5) 

The conformal metric is assumed to be regular at 
in compactificd coordinates. This implies that $1 = 
at since, in compactificd coordinates, the physical 
metric is singular there. Comparing Eq. ([5]) with the 
more widely used definition g^ = "^^gij , we see that f2 = 
-0^^. The advantage of the definition (O is that is 
finite at J"^ . The equations in the rest of this section 
mirror the standard conformal method [ssl . [39j , with the 
replacement f2 <-> ^/;"". 

The trace-free extrinsic curvature is defined as 

= i^y - i g^^K, (6) 
3 

where A"'^ = g'^^g^'^Kira- The physical trace- free extrin- 
sic curvature A^^ and its conformal counterpart A^^ are 
related by 

^y=17iy, A^ ^Vc" AK (7) 

Substituting Eq. ^ into Eq. (g]) and using the CMC 
condition VjiC = gives 

Vj^^^' = 0. (8) 
This can be rewritten as 

Vc'ij^A'^ = 0, (9) 

where Vi is the covariant derivative with respect to ijij. 
Because of the CMC condition, the momentum con- 
straint has decoupled from the Hamiltonian constraint, 
and can be solved first. Standard methods for solving 
the momentum constraint are given in Ref. [40| . In this 
paper, we assume a fiat conformal metric gtj. 

Substituting the solution A*-' of the momentum con- 
straint into the Hamiltonian constraint ([5]) and express- 
ing it in terms of conformal quantities, one finds 
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Here denotes the covariant Laplacian with respect 
to gij and TZ is the scalar curvature of Py . Note that 
some terms in Eq. (jlOp arc singular at where 51 = 
0. Assuming Aij is finite at .y'^, any regular solution 
satisfying the boundary condition 



0; it intersects past null infinity. We arc interested in 
hypersurfaces approaching future null-infinity, and will 
therefore require 



K>0 



(16) 







must also satisfy 



(11) 



(12) 



throughout this paper. 

The constant C determines whether the Et hypersur- 
face intersects the black hole horizon or the white hole 
horizon. Considering radial null rays on the horizon, 
r — 2M. we find that for 



B. The Schwarzschild black hole in CMC slicing St intersects the black hole horizon (i.e. 



(4), 



(17) 



< for 



Let us now discuss this initial value formalism in the 
particular case of a Schwarzschild black hole, recasting 
already established results [U, \^ in our language. 
In spherical coordinates (r, 6, 4>, t) , with t constant on 
CMC slices and r the areal radius, the metric is given 
by 0,113, m 



1 



2M 



1 



2a 



dr + -Kdr^ -dtdr 



(13) 



sin^ e c 



where the functions / ~ f{r) and a = a{r) are 



/(r)= (1-^ + a 



1/2 



a{r) = - 2- (14) 



With these definitions, a(r) — '^'n'", where ^'^'n'' is the 
radial component of the hypersurface-normal '■^^n^. The 
constant M is the mass of the black hole. The constant 
C represents an additional one-parameter degree of free- 
dom in the choice of spherically symmetric CMC hyper- 
surfaces in the Schwarzschild metric. The lapse-function 
a of the metric (fT5|) equals /, and the only non- vanishing 
component of the shift is = —af. 

The coordinate transformation between Schwarzschild 
coordinates (r, 9, (j), T) and the CMC coordinates 
{r,9,(f>,t) (again, see HI S El ) is 



t 



T 



a{r) 



2M^ 



fir) 



-dr 



(15) 



dr 



1 



2M^ 



a{r) 



1 



2M^ 



fir) 



-dr. 



where u is the Eddington-Finkelstcin retarded null co- 
ordinate. With K > 0, constants of integration can be 
chosen so that t — > u as r — > oo. 

The hypersurfaces S( of constant coordinate value t 
will play a central role in this paper. For K ^ 0, these 
hypersurfaces extend to space-like infinity. For K 0,'Et 
becomes asymptotically null for large radius. If K > 0, 
this hypersurface intersects future null-infinity, because 
1// ^ and a/f^ +1 in the limit r oo. For K < 



r < 2AI). If the inequality is reversed, the hypersur- 
face enters the white hole region, where the light cone is 
tilted toward increasing r, and excision is is not possible 
without allowing causal propagation from the excision 
boundary to the interior of the computational domain. 

As with any spherically symmetric metric, a radial co- 
ordinate transformation r ^ R ~ Rir) can be used to 
make the spatial sector of Eq. p3p conformally flat, i.e. 



Here, 



r 



(18) 



and R is determined by the ordinary differential equation 

dR _ R 
dr rf 

oo. 



(19) 



Because f ^ r for large r, R remains finite as r 
Denoting its limiting value by we find 



Rir) 



exp 



dr' 
r'fir') 



(20) 



As r — s- oo (or equivalently R R+), fl ^ 0, in agree- 
ment with the boundary condition Eq. (|lip . 

Finally, transforming (R, 6, 0) to Cartesian coordinates 
x\ we can express the space-time metric Eq. (fT3)) as 



ds^ = n- 



a^'dt^ + S,j{dx' + f3'dt){dx^+(3^dt) . (21) 



Here the conformal lapse a and conformal shift are 
given by 

a^flf and ff^-flan\ (22) 

where = x^/R. Because il ^ R+/r whereas f ^ a ^ r 
as r ^ oo, a and /?' are finite at J^+. Therefore, the 
conformal space-time metric inside the square brackets 
in Eq. (PT|) is regular, in addition to its spatial part being 
flat. 

The trace-free extrinsic curvature of the hypersur- 
face in the coordinates (t, a;*) takes the form A^^ = fl^A^^ , 
with 



i?3 



37i'n^ - 6 



(23) 
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The sign difference in Eq. (|23p relative to earlier work 
(e.g. Eq. (52d) of [111) arises because of the different sign- 
convention for Aij [see the discussion after Eq. ([2])] . The 
conformal trace-free extrinsic curvature satisfies the mo- 
mentum constraint Eq. ([9]). Indices on conformal spatial 
tensors such as A^^ are raised or lowered with the con- 
formal spatial metric; in our coordinates, this is simply 
the Kronecker delta, so the components of such tensors 
are identical irrespective of the index-location. 

We finish this section by noting that Eq. ([20]) deter- 
mines only i?(r)/i?+; the freedom remains to rescale R 
by a real constant 77, 



This rcscaling induces further rcscalings. 



n 

An 



7]-^A'', 
i]-'^Aij. 



(24a) 



(24b) 
(24c) 
(24d) 



Eqs. represent a coordinate transformation and do 
not affect the physical initial data. For the zero mean 
curvature case, K = 0, the hypersurface asymptotes to 
spatial infinity, where it is natural to impose the bound- 
ary condition 17— >1, soi?/r^lasi?— >oo. With 
K > 0, the coordinate scale is set by the arbitrary choice 
of the value of R at future null infinity, R^. Only af- 
ter the Hamiltonian constraint is solved can R and be 
rescaled to make the maximum value of equal one. If 
K is very close to zero, this rescaling makes H very close 
to the K = solution almost everywhere. 



C. Minimal Surfaces & Trumpets 

Later in this paper, we will construct black hole initial 
data with minimal surface boundary conditions in the 
interior of the black hole(s). As preparation, let us dis- 
cuss the presence and location of minimal surfaces in the 
hypersurfaces of the metric Eq. (fT3|) . Our treatment 
here extends earlier similar discussions [4l|, [Z^, . 

The sphere r =constant is a minimal 2-sphere within 

if the function / defined in Eq. (fT4)) vanishes. This can 
be seen most easily by noting that r is the arcal radius, 
and that from Eq. (|19p . dr/dR = fr/R, which vanishes 
for / = 0. For any values M > 0, K and C, /(r) > 
for r > 2M. Furthermore, the inequalities Eqs. and 
([T7]) imply that / is strictly positive, / > 0, at the black 
hole horizon r ~ rn = 2M. Therefore, a minimal surface 
in St will always lie inside the horizon, if it exists. 

Vacuum general relativity possesses a rescaling free- 
dom, and this freedom will be inherited by the func- 
tion / and any equations that describe minimal surfaces. 
The most common way to incorporate this freedom is to 
rescale all dimensionful quantities by the mass M. In 
particular, Eq. (|14p can be written in terms of the di- 




FIG. 1: Radial function /(r) for different values of K and C. 
Plotted is the square for KAI = 1/10 and several different 
choices for C. Each curve is labeled by its value of C/M^ and 
its value for K^C. The bottom axis shows radius in units of 
r/M, and the top axis in units of Kr. 



mensionless variables r/M, KM and C /M'^ as 



f = l- 



7M 



{KM){r/M) C/AI^ 



(r/A/)2 



(25) 



Using M to make dimensionless variables is fine for 
Schwarzschild, but in the absence of spherical symme- 
try, the mass M is not known in advance. On CMC 
hypersurfaces, K is a free parameter and can be used 
to form an alternative and more widely applicable set of 
dimensionless variables, Kr, K^\/C, and KM, such that 



2KM 

/ = 1 + 

Kr 

A minimal surface at radius 



Kr 
~3~ 



K'C 



{Krf 



2KM 



Kr„ 



satisfies 



K^C 

{Kr^^sY 



= 0. 



(26) 



(27) 



This equation will play a central role in the remainder of 
this section. 

Fig. [1] plots p for different values of C/M"^ for a fixed 
value of KM. For large C/M^ (or equivalently K'^C), 
no root and no minimal surface exists. At some critical 
value Ct/M^, f touches zero, indicated by the filled cir- 
cle. For C/M^ < Ct/A'/|, minimal surfaces exist at radii 
Tms/M, as indicated by the open circles in Fig. [TJ The 
critical value Ct/M^ delineates the region of parameters 
for which contains a minimal surface. At this critical 
point, / varies linearly in r — rx, passing through zero at 
rx. From Eq. (|13p . the radial proper separation within 
the shce is ds = dr/f. Because / vanishes linearly at rx, 
this point is an infinite proper distance away from any 
point r > rj : this configuration is often called a trum- 
pet [H, [2^, [s^, [3l|. In contrast, away from the critical 
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point (i.e. at the open circles in Fig. [T]), / approaches 
zero proportionally to — r„is, and the proper distance 
to the minimal surface is finite. 

For trumpets, = and 9r.f ^ = 0; the parameter val- 
ues that lead to trumpets can be written down in terms 
of the arcal radius of the trumpet, rx: 



KtMt 



2rT/MT - 3 



(rx /Mt ) y/{rT /Mt) (2 - /Mt) ' 

(rT/MT)2(3- tt/Mt) 
3V(rT/A/T)(2-rT/MT)' 



(28a) 
(28b) 



where 3/2 < tt/Mt < 2. 

The preceding discussions determine the region of pa- 
rameters M, iC, C, for which minimal surfaces exist. Tak- 
ing the scaling invariance into account, this region can be 
represented on a two-dimensional plot as given in Fig. [2l 
The blue solid and red dashed lines in the top panel, for 
instance, are given by a parametric plot of Eqs. (|28p and 
Eq. PT|) . The unshaded wedge-shaped region between 
these two lines represents the allowed parameter choices 
which lead to a CMC hypersurface containing a minimal 
surface and intersecting the black hole. 

Alternatively, one can compute trumpet-configurations 
using the dimensionless variables indicated in Eq. (|26p . 
From the equations p = and drp = 0, one elimi- 
nates KM to obtain a third order polynomial that re- 

1/2. 

lates K^rT and K^^Cr^ . This polynomial has only one 
positive real root, which for K^Ct < 2/3 can be written 
as 



cos [i arccos (|7^|,Ct) 



(29) 



For K^Ct > 2/3, the trigonometric functions in Eq. ((29|) 
should be replaced by their hyperbolic counterparts. To 
find the trumpet solution, substitute K^Tt back into 
Eq. ([261): 



KtMt 



K^rT 



(30) 



Substituting Kttt from Eq. ([29]) . KtMt is a function of 
K^Ct alone. Finally, dividing Eq. ^ by Eq. ^ yields 
the value of vt/Mt- All these parameters for trumpet 
hypersurfaces are plotted in Fig. \3\ To make easy con- 
tact with dimensionless quantities normalized by M , the 
top horizontal axis of this plot is labeled by KM. For 
KC^^^ <C 1, the data plotted in the lower panel of Fig. [3] 

1 /2 

is proportional to XtC't ■ 



Kttt - Z^'^KtCt'^ + O ((/V 2 ^t)'/') , (31) 
KtMt = -^KtC^^ + O [K^Ct) . (32) 



For given choices K and K^C, minimal surfaces only 
exist at radii r^s larger than tt given by Eq. (|^^ . This 
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hypersurface 
intersects white hole 
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1 1 i-i - ■ ■ - I- 
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FIG. 2: Parameter choices that result in hypersurfaces Et con- 
taining a minimal surface. The two panels correspond to two 
different choices of dimensionless variables. Parameter values 
on the blue solid lines represent trumpet configurations. 




0.001 0.01 



FIG. 3: Properties of the trumpet hypersurfaces, 
parametrized by the dimensionless parameter KC^^^ . 
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iTrms shows that this hnc is parametrized by 



0.0001 0.001 0.01 0.1 1 

FIG. 4: Properties of CMC hypersurfaces viewed in the 
Rms/R+ vs. KC^^^ plane. CMC hypersurfaces with a min- 
imal surface must lie below the thick black line, and as this 
line is approached, the minimal surface approaches the black 
hole horizon. The dashed red lines are lines of constant MK, 
with values given by the red numbers next to the lines. The 
thin blue lines represent constant values of 7 = rms/rn- The 
shaded areas are contours of constant value of Rms/RAW, the 
shade of grey changes, from top to bottom, at values 0.9, 0.5, 
0.1, 0.01, and 0.001. Trumpet hypersurfaces represent the 
limit Rms/R+ 0. 



is indicated by the solid blue line in the lower panel of 
Fig.H 

To obtain an upper limit on Kr^g, we recall that all 
minimal surfaces must lie inside the horizon, r^g < fH = 
2M. Combining this with Eq. pT|) results in 



Kr n,s < 2MK < {3K^C) 



1/3 



(33) 



which is indicated by the red dashed line in the lower 
panel of Fig. [2] 

Finally, Fig. 2] presents another view of the two- 
dimensional set of "good" parameter choices that we first 
indicated in Fig. [2] As in the lower panel of Fig. [2l we 
shall use KC^^^ to parametrize the horizontal axis. How- 
ever, we use Eq. (|20p in the form 



In 



l/{Kr„ 



du 



■ - 2KMu^ + - K^Cu^y 

(34) 

to convert the vertical axis to the ratio Rnw/R+ of con- 
formal radius of the minimal surface and the conformal 
radius of J^+. Fig. [4] shows the following data: First, the 
thick black line corresponds to parameters for which 
the minimal surface coincides with the black hole hori- 
zon. Setting Kr^is — 2KM in Eq. (P7|) and solving for 



Kr„ 



(3A'2C)i/3, 



(35) 



with /vrms mapped to i?,„s/i?+ by Eq. ([M)) . The red 
dashed lines are contours of constant values KM . Each 
of these lines is obtained by keeping KM fixed and vary- 
ing Kr^is between its lower bound, the trumpet value 
/^T'^T (obtained from inverting Eq. (|28ap ) and its max- 
imal value 2KM. For each choice of KM and Kr, 
Eq. ([?7|) is solved for K^C, and the resulting data plot- 
ted as a parametric plot. Finally, the thin blue lines in 
Fig. m represent lines of constant ratio Kr^s/ {2KM) = 
fms/fH = 7, he. lines where the areal radius of the min- 
imal surface is a fixed fraction of the areal radius of the 
horizon. Replacing Kr^s by 2'yKM in Eq. (P?]) . we can 
solve for KC^/'^ as a function of KM. The thin blue lines 
are then obtained as a parametric plot {KC^^^ , Kr) as 
KM is varied. 

Trumpet initial conditions are obtained from Fig. [3] 
through the limit R,^s/R+ 0, i.e. by going "down". 
Note that the red KM =constant contours become ver- 
tical in this limit. Their value as a function of KC^^^ is 
then given in the lower panel of Fig. O 

The significance of the axes employed in Fig. 0] is that 
both axes represent quantities that are freely specifiable 
when constructing CMC hypersurfaces within the initial 
value formalism of general relativity. One use of Fig. 2] 
is to first pick values of K and M, fixing a particular 
contour of KM. One then chooses how trumpet-like the 
initial conditions should be; that is, how far to go down 
along the contour. Alternatively, one can choose a certain 
ratio 7 = Ti-asl^H- In either case, one can then read off 
the corresponding values of Rms/R+ and K^/C to get 
the remaining initial value parameters. 



D. CMC initial data for single black holes 

The properties of CMC slices of the Schwarzschild 
spacetime, as described in Sees. Ill Bl and III C\ mesh 
nicely with the conformal method of solving the Einstein 
constraint equations, which was outlined in Sec. Ill Al 
We shall first discuss the spherically symmetric case, 
for which there is a one-to-one correspondence. Subse- 
quently, we will generalize to a single spinning or boosted 
black hole. 



1. Spherical symmetry 

The CMC metric is conformally flat, so we shall use a 
flat conformal metric for the initial value problem. 



9ij — Ik 



(36) 



where /y is the flat space metric. The extrinsic curvature 
of the CMC metric has the correct scaling with confor- 
mal factor [compare Eqs. ([7]) and so we shall adopt 
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Eq. p3|) as the freely specifiable tracefree extrinsic cur- 
vature, with the constant C yet to be determined. The 
radial coordinate R ranges from a finite value repre- 
senting to some smaller value, for instance i?ins at 
a minimal surface (assuming a minimal surface exists), 
so we shall adopt a computational domain with inner 
radius i?i and outer radius i?2- At J^'^ , the conformal 
factor vanishes, resulting in the boundary condition 



n = 0, 



(37) 



which identifies i?2 with At the inner boundary, we 
shall impose a minimal surface boundary condition, 



dR 



n 



(38) 



so that Ri will coincide with Rms- 

With the choices Eq. ((36)) - p8l) . we are now left with 
choosing the four numbers {A', C, i?2}- Solution of 
the Hamiltonian constraint Eq. ((T0|) will then result 
in a complete initial data set with a certain mass M. 
Fig. |4] is useful for informed choices for the numbers 
{K,C, Ri, R2}- For instance, we can first decide on a 
mass M (say, M = 1) and a mean curvature K (say, 
K = 0.01). This selects one of the red-dashed curves in 
Fig. m We can now choose a suitable value of r„^s/rH 
by considering the intersection of the red dashed lines 
with the blue contours (say, r,„s/?'H = 0.9), and read 
off the values for KC^^^ and i?nis/-R+ (in our example 
KC^/^ = 0.0105, Rn,s/R+ = 0.000554), which determine 
the values for C and i?i/i?2. An overall scaling of Ri 
and i?2 remains, because the coordinate transformation 
r ~* R for CMC slices is determined only up to an overall 
rcscaling (see the discussion after Eq. (PO)) ). Thus, we are 
free to set, for instance, Ri — \. 



2. Single black hole with spin & boost 

A relatively simple class of non-spherically symmetric 
initial data on maximal slices, with K = was pro- 
posed by Bowen and York [l^. El]- It assumes confor- 
mal flatness of the spatial metric and a solution Aij of 
the conformal momentum constraint Eq. ^ character- 
ized by a "spin" vector S** and two "boost" vectors P*, 
Q*. On the asymptotically flat maximal slices, with the 
boundary condition ^ 1 at spatial infinity, S*' is in fact 
the physical angular momentum and P' is the physical 
linear momentum of the system, as defined by Arnowitt- 
Deser-Misncr (ADM) surface integrals at spatial infin- 
ity [l^]. The second boost vector Q* is introduced to 
allow for inversion symmetry about a minimal surface, 
and can be thought of as the three-momentum of the 
black hole as viewed from the asymptotically flat space 
on the other side of the Einstein-Rosen bridge associated 
with the minimal surface. 

Since Eq. ^ is linear and identical on maximal 
and CMC slices, we can add the Bowen- York terms to 



Eq. ((231), with the result^: 
C 



(39) 



2P2 



Here is a unit three-vector in the outward radial di- 
rection. The coefficient C of the spherically symmetric 
first term has been normally taken to be zero in papers 
on the initial value problem on maximal hypersurfaces. 
In this paper, we refer to Aij given above in Eq. ([M)) as 
the generalized Bowen-York solution. 

Note that A*-' is not invariant under the rescaling of R 
(see Eqs. ([24])). A scale-invariant effective source term in 
Eq. Uni) is 



R'^A.jA'^ = W{R,9,'fi). 



(40) 



The properties of W are important for determining ques- 
tions such as the inversion symmetry of the hypersurface 
about the minimal surface (see the Appendix). Also, the 
invariance of W under the rescaling freedom R i]R 
implies that the parameters must scale as 

C->C, S'^S\ P'-*r)-^P\ Q'^f]Q\ (41) 

The physical interpretation of the Bowen-York param- 
eters is not necessarily the same on CMC hypersurfaces 
as on maximal hypersurfaces. CMC hypersurfaces are 
not asymptotically flat. Identification of the physical en- 
ergy, linear momentum, and angular momentum of the 
system on asymptotically null hypersurfaces is a non- 
trivial matter, in general [45j . In particular, the scaling 
dependence of the boost parameters means that these 
cannot be interpreted as physical momenta. 

One can argue that in the limit of small K and C 
{KM < 1 and C - {8M'^/3){KM) < 1), the geom- 
etry in the vicinity of the black hole should be simi- 
lar to that of solutions of the zero mean curvature ini- 
tial value problem found, for example, by Cook 0, HJ. 
For hyperboloidal slices, the conformal factor is gen- 
erally approximately constant at intermediate distances, 
flM < i? < 0/A'. (As can be seen in Fig. [71 is smaller 
near and inside the black hole, where fl R/M, and de- 
creases toward zero approaching ^l < KR.) These 
intermediate distances are sufficiently close to the black 
hole that the CMC slice still resembles a maximal slice, 
but far enough away to be considered in the asymptotic 
regime. If R is rescaled such that sa 1 in this interme- 
diate regime, then the ADM formulas for energy, linear 



^ The sign differences between Eq. H39| l and earlier papers arise 
because of the different sign convention for Aij (see Eq. ((2]l ) . 
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momentum, and angular momentum should be at least 
roughly valid. Therefore, one might identify the scaling 
invariant Jlmax-P' as a quasi-local linear momentum and 
S"* as a quasi-local angular momentum. However, these 
"quasi-local" values may not match the correct physical 
values at future null infinity if gravitational radiation is 
present outside the plateau region. As a surrogate for 
the mass of the black hole, we use the "irreducible mass" 
Mirr, defined in the usual way from the area of the ap- 
parent horizon Aahi 

If the initial data is axisymmetric, then the physical an- 
gular momentum can be calculated precisely using stan- 
dard techniques. The generalized Bowen-York solution 
for a single black hole is axisymmetric provided that the 
boost is zero or aligned with the spin vector. Then the co- 
ordinates can be chosen so that only the z components of 
the spin and boost vectors are non-zero. The solution for 
the conformal factor will be axisymmetric for our mini- 
mal surface inner boundary condition, since the minimal 
surface is assumed to be a coordinate sphere. Choose 
spherical polar coordinates (i?, 0, (p) in the conformal 
flat space, so that the axial Killing vector ^ C = d/dip. 
Then the Komar angular momentum within a coordinate 
sphere is 

I I aVg^'^'^dedp (43) 

J^TI" Jo Jo 

1 r'2'K 

In much of the earlier work on maximal hypersur- 
faces, considerable emphasis has been placed on ob- 
taining inversion-symmetric solutions of the initial value 
problem (see [13, 0, \^). Given a minimal surface at 
R = Rms, inversion symmetry requires that Q{R, 9, ip) ^ 
{R^^/Rfn{RlJR,e,ip). Solutions of Eq. ^ for n 
with minimal surface boundary conditions can be contin- 
ued with inversion symmetry to i? < R-ms if and only if W 
as defined in Eq. (140]) is inversion symmetric. For C ~ 0, 
the usual story is that the boost terms are inversion- 
symmetric if = ±i?^jgP*, with only the minus sign 
applicable if spin is also present. We show in the Ap- 
pendix that with C ^ 0, inversion symmetry requires 
the plus sign when only boost is present, and that no in- 
version symmetry is possible with both boost and spin 
unless the boost and spin vectors are co-linear. Inver- 
sion symmetry is desirable primarily for simplifying ex- 
cision boundary conditions during evolution. On CMC 
hypersurfaces, the generic absence of inversion symmetry 
requires rethinking how excision will be handled, and it 
may be desirable to just set = 0, as also done in [23t . 

The location of the apparent horizon relative to the 
minimal surface is an important issue. The inner bound- 
ary of the computational domain should be inside the 



apparent horizon. In spherical symmetry, the apparent 
horizon will always be outside or on the minimal surface, 
provided that a solution of the Hamiltonian constraint 
exists. However, in the presence of a boost, the appar- 
ent horizon can straddle the minimal surface [47| . Care 
should be taken in the choice of R^s as the boost is in- 
creased, in order to avoid this problem. 



E. Multi-black hole solutions 

The goal of this section is to construct a CMC hyper- 
surface with N black holes, with masses (approximately) 
A/q,, Bowen-York parameters Pa and Sa at coordinate 
locations c^. Here a labels the black holes. 

The setup for single black holes in Sec. Ill Dl can be 
generalized by choosing one excision boundary for each 
black hole, with radius Ra centered at Ca- Because the 
momentum constraint is linear, the extrinsic curvature 
can be taken as the superposition of N copies of Eq. (|39p. 
each one centered at the appropriate Cq,: 

= E (^S"" + ^S"^ + + ^%') (44) 

a 

If the black holes are sufficiently widely separated, and 
if the outer boundary is sufficiently far away, we expect 
that close to each of these black holes, the solution is a 
perturbation of the single black hole case. 

In the asymptotically flat case, the conformal factor 
is close to unity, except very close to each black hole. 
Therefore, the coordinate distance |cq — cp\ between the 
black holes a and /3 is a convenient and reasonably ac- 
curate approximation of the proper separation between 
the black hole horizons. Because of the rescaling freedom 
discussed in Eqs. ([M)) . may not be close to unity on 
the hyperboloidal slices considered here, and therefore 
the coordinate distance may deviate signiflcantly from 
the proper separation. 

There is only one global value K and one value i?+ for 
the whole multi black-hole configuration, whereas each 
black hole has its "own" constants Ca , Ra , as well as P^ 
and S'a- Therefore, the interesting question arises of how 
to use Ca and Ra to control properties of the individual 
black holes, for instance their masses, given fixed values 
for K and R+. Assuming that the presence of and 
Sa will only mildly perturb the case of the spherically 
symmetric black hole, we can use Fig. |4] to address this 
question. A given K and a (desired) value for M = 
Ma places the solution on a particular KM = constant 
contour. Given a desired value for 7 = r^s/i^H, a unique 
point in this figure is determined, and one can read off 
Rms/R+ and KC^^'^ and then compute Ra = Rms and 

Ca = C. 

This procedure can be simplified in the particularly 
interesting limit KM <C 1. Consider a minimal surface 
with 3/4 < Vms/fH = 7 < 1, and with a given value 
of KM. Substituting KM and Kr^s = 2jKM into 
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0.08 



0.04 



FIG. 5: The functions gc and gn, which are the asymptotic 
values of C/M^ and i?ms/(-R+ KM) in the hmit KM 0. 
These functions are useful for choosing excision radius and 
C when constructing CMC initial data for black holes with 
given K, M and outer boundary radius R+. 



blows up at i? = 0. This is a very reasonable condition 
which follows automatically from inversion symmetry in 
the limit i?nis (see the Appendix), and reflects the 
fact that the other side of the Einstein-Rosen bridge is 
infinitely far away from any point with R > 0. 

We begin by rewriting Eq. (|10[) with U = Vl/R as the 
dependent variable. The new form of the equation is 



R 



.dHJ_ 



4i? 



3 

2U 



U + R 



dR 
dU_ 
dR 



2U + AU 



1 



-WW 



(47) 



where A is the Laplacian operator and V the gradient 
operator on the unit two-sphere, and W is defined by 
Eq. (|40|l . We have assumed a conformally flat spatial 
metric. 

Now let U = Uo{e,Lp) + R'^Ui{e,ip) + ... aiidW = 
Wo(0, tp) + RWi{6, (f) + . . .. From the expression for W 
in Eq. (jA.2p . we sec that 



Wo = 6(C2-h3sin2(V;) S'S^) 



(48) 



Eq. (P7)) . one can solve for A'^C and compute C. The 
result is 



KM < 1 



(45) 



with gc{l) — 47^7 — 72. We also find numerically that 



QRi-f) KM, 



KM < 1. 



(46) 



Therefore, if one knows the coefficients gc and g^, for the 
desired ratio 7 — r-^s/ru, one can immediately compute 
the values for C and Rnis/R+ from Eqs. and P5|) . 
These coefficients gc{l) and 5^(7) arc plotted in Fig. 



F. Trumpet Inner Boundary 

In the absence of spherical symmetry, a trumpet solu- 
tion of the Hamiltonian constraint equation can be un- 
derstood as the limit i?,us — *■ 0. That is, the Hamiltonian 
constraint Eq. (fTU|) is solved with the boundary condition 
dfl/dR = n/ R a.t R = 0. The consequence of this bound- 
ary condition is that = and U/R is finite at i? = 0, 
which in turn means that the proper distance from finite 
i? to i? = 0, dR/Q,, is infinite. In this section, we 
show how the singularities in the equation determine the 
non-singular solution at the trumpet inner boundary and 
derive the behavior of the solution close to the trumpet 
boundary. 

Note that a necessary condition for a non-singular so- 
lution is that the second "boost" vector = 0. Other- 
wise, the W source term in the Hamiltonian constraint 



for any boost P', where tjj is the angle with the spin 
direction, equal to 9 if the spin is along the polar axis. 
Unless the boost is non-zero, W = Wq at all R and Wi = 
0. 

In zeroth order, we get 



1 



3 
2C^ 



VC/o • VC/o 



WoU^ 



(49) 

which has a unique solution regular everywhere on the 
unit sphere for any A' > 0, any value of C > 0, and 
any spin vector S*'. Uniqueness can be demonstrated 
using a method of Moncrief [s^l applied to the quasi- 
linear form of the equation obtained by the change of 
variable Uq ^l/V^. 

In the absence of spin, Uq is independent of angle and 
Uq is the solution of the cubic equation 



The only positive real root if 3Ar^C/2 < 1 is 



(50) 



C/o- = ^cos 



-K^C 



(51) 



The trigonometric functions are replaced by the corre- 
sponding hyperbolic functions if 3Ar^C/2 > 1. 

The ncxt-to-lcading terms in Eq. (j47|) give an equation 
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for Ui: 



AUi 



— ur- 



1 - ^WoU^ 



(52) 



Uo 



If Wi 7^ 0, the solution of the inhomogeneous equation 
with a = 1 gives the leading contribution to Ui . There is 
also a unique lowest value of a for which the homogeneous 
equation has a non-trivial solution regular everywhere 
on the unit sphere. This solution to the homogeneous 
equation, times R", will, with a coefficient undetermined 
by the trumpet boundary condition, contribute to U—Uq. 
The coefficient is fixed by the requirement that the global 
solution for fi satisfy the $1 = boundary condition at 
future null infinity. 

If the spin is zero, Wq ~ 6C^ and the homogeneous a 
is the solution of the algebraic equation 



a 



1 = 9C^L/o = 4cos2 



1 



K^C 



(53) 



In the range < 3K'^C/2 < 1 of most interest, Eq. ^ 
implies -^2 < a < ^/3. For larger K^C, the trigonomet- 
ric functions are replaced by hyperbolic functions and a 
continues to increase. 

In practical terms, there is very little difference be- 
tween a solution satisfying the exact trumpet bound- 
ary condition and a solution satisfying the minimal sur- 
face boundary condition with a very small, but non-zero, 
Rnis- Very small means that i?i„s/i?AH <C 1. For the 
Schwarzschild case, i?,„s/i?+ should be far below the 
heavy black line in Fig. [H 

Finally, all of this discussion has been in the context of 
single black holes. With multiple black holes, each trum- 
pet boundary must be treated separately and matched 
to the global solution on a surface surrounding the black 
hole. The analysis right at the trumpet boundary is not 
affected by the presence of other black holes, since the 
R^ factor in W kills the finite contribution of the other 
black holes to the conformal traceless extrinsic curvature 
at i? = 0. 



III. NUMERICAL RESULTS 

In this section, we numerically construct a variety 
of hyperboloidal initial data sets using the generalized 
Bowen-York solution. These results are obtained with a 
pseudo-spectral elliptic solver that is part of the Spectral 
Einstein Code, SpEC. This solver is described in detail in 
Ref. [131 ■ The desired solution is expanded in terms of 
spherical harmonics and Chcbyshcv polynomials. Trun- 
cation at some finite expansion order results in an alge- 
braic system of equations for the expansion coefficients 



or, equivalcntly, for the values of the solution at the col- 
location points. This system is solved with a Newton- 
Raphson technique, employing the preconditioned gener- 
alized minimal residual method (GMRES) [5l| to solve 
the linearized system of equations at each iteration us- 
ing the software package PETSc [H,!!!,!!!. The SpEC 
elliptic solver has been used on a wide variety of formu- 
lations of the initial value problem (see, e.g. [H, |5^]), 
including puncture initial data [13, [ll] (which also uses 
the Bowen-York extrinsic curvature). Below, we present 
results of convergence tests of our initial data. 



A. Spherical Symmetry 

As a first test, wc reproduce the analytically known 
spherically symmetric solutions discussed in Sec. Ill Bl us- 
ing the numerical approach described in Sec. HID II We 
choose M = 0.85, K = 0.1, r^s/rH = 0.8 and R+ = 100. 
The relations shown in Fig. |4] then imply C = 1.0086 
and _Rnis = 0.127. From these numbers, only K, C, 
Rms ~ Ri, and = i?2 are used in the numerical 
solution; the other numbers arc used only when com- 
puting the analytical solution with which to compare. 
Fig. [5] shows convergence of the numerical solution to 
the analytic solution. The solid lines plot the differences 
between the numerically determined and the analytic 
solution of Eq. ([34]) computed with Mathematica. As we 
increase the resolution of the elliptic solver, we find ex- 
ponential convergence to the analytic solution. For the 
generic examples considered later in this paper (which 
include spin, boost, and two black holes), no analytic so- 
lutions are known. Therefore, in Fig. [SI we also present 
an estimate of the numerical error which does not rely 
on knowledge of the analytic solution. Specifically, the 
dotted lines show the differences between the numeric 
solutions at two successive resolutions. As can be seen, 
these track very closely the error obtained from compar- 
ing the lower resolution run to the analytic solution. 

In our second example, we explore a solution which is 
very close to the trumpet configuration. Recall that for 
a trumpet, for R close to i? = 0, the solution behaves as 
r2 = UqR with Uo given in Eq. ([5T|) . We choose parame- 
ters C = 1, K = 0.1, = 10"^, and R+ = 100 which, 
as can be seen from Fig. [4l result in an inner boundary 
which is very close to the trumpet limit. Because applica- 
tion of the minimal surface condition in this case proved 
numerically problematic (presumably due to dividing by 
the very small number i?ms), we use the Dirichlet condi- 
tion = UoRyns at the inner boundary. 

The numerical solution of the Hamiltonian constraint 
equation for this example is shown in Fig. [7| One sees 
that in the region 10^^ < ^ < 1, the conformal factor 
is proportional to R, with Uq the constant of proportion- 
ality. Furthermore, within this range of conformal ra- 
dius, the proper area of coordinate spheres (47rr^, where 
r is the Schwarzschild radius) is approximately constant, 
which is consistent with the long cylinder of the trumpet. 
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FIG. 6: Convergence of the numerical solution of the Hamil- 
tonian constraint for a Schwarzschild black hole. Plotted are 
five different resolutions Nr, where Nr is the number of ra- 
dial collocation points. Solid lines show the difference from 
the analytic solution, dotted lines the difference from the nu- 
merical solution at the next higher resolution. The highest 
resolution has Nr = 104. 
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FIG. 7: Numerical solution for a Schwarzschild black hole 
with the inner boundary very close to the trumpet limit. The 
dashed red line shows the conformal factor f2, and the solid 
blue line shows the Schwarzschild radius (calculated from the 
proper area of coordinate spheres) as a function of the con- 
formal radius. The dotted black line is UqK, which equals Q. 
for a trumpet, with Uq given in Eq. (|51|l. The vertical dashed 
line locates the apparent horizon. 



The location of the apparent horizon rjj = 2M = 1.69 is 
shown vertical dashed Unc in this figure. 



B. Single spinning black hole 

Here, we construct a single spinning black hole, with 
no boosts. We take the Bowen-York spin parameter 
S"* = (0, 0, S) and solve the Hamiltonian constraint for 
the conformal factor O, with varying S. The solution is 
axisymmctric, so S"* represents the total angular momcn- 
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FIG. 8: Irreducible mass versus dimensionless spin for a single 
black hole spinning around the z-axis. 



tum of the black hole (see Sec. HID 2p . 

In the absence of boosts, W (Eq. (001)) reduces to Wq 
as defined by Eq. with ij^ = 9. The radial behavior 
of f2 will be rather similar to the spherically symmetric 
solution with the same parameters as long as the spheri- 
cally symmetric is replaced by the solid angle average 
of VF/6, which we denote by 



25' 



(54) 



We use S/CcS as a dimensionless measure of the impor- 
tance of spin. Note that < S/Cqs < 1/^2 as S/C 
varies from zero to infinity. We find that the irreducible 
mass varies less with spin keeping Ccs constant than 
when keeping C constant, particularly for large spins, 
as shown in Fig. [5] The constant values of C and CcS 
arc 1.0086, with K = 0.1, R,^s = 0.127, and R+ = 100. 

In Fig. [9l wc study how the non-zero spin distorts the 
intrinsic geometry of the apparent horizon and compare 
the results to an analogous distortion computed from 
the analytic Kerr solution. The solid lines of Fig. [5] 
show the maximum and minimum of the Ricci scalar 
^"^^TZMf^^ computed from the 2-mctric induced on the ap- 
parent horizon of the spinning black hole. The dashed 
lines show the maximum and minimum values calculated 
from the analytic Kerr solution, taken from Eq. (Bl) of 
Ref. [111. Note that deviations from 0.5 are deviations 
from a spherical geometry. We see that the apparent 
horizon distortion is much less for our conformally flat 
initial data than it is for Kerr. The CMC data plotted in 
Fig. [5] is the same shown in Fig. [51 with the CMC curves 
terminating at max(S'/Coff) = 1/^/2. 

The horizontal axis in Fig. [5] is the spin-cxtrcmality 
parameter 



c 



(55) 
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FIG. 9: Dimensionless scalar curvature of the apparent hori- 
zon versus spin for a single black hole. The solid lines are the 
maximum and minimum numerical values computed on CMC 
slices. The dashed lines are the maocimum and minimum an- 
alytic values for a Kerr black hole (these lines continue to the 
maximal Kerr value C = 1: where the maximum and minimum 
are 2 and —1/2, respectively). 



as introduced in Ref. [58j . A maximally spinning Kerr 
blaek hole has C = Ij and the CMC-sequence considered 
here allows values as large as C ~ 0.78. In Sec. |IVl we 
shall place this number into the context of results on zero 
mean curvature slices. 



C. Single boosted black hole 

Next, we construct single, non-spinning, boosted black 
holes. We shall vary and shall choose = +l(^^^P'^ , 
in order to make the black hole spacetime inversion sym- 
metric (see the Appendix for details). As we vary the 
boost, we keep the irreducible mass of the constructed 
black holes constant by a suitable choice of i?ms- Specif- 
ically, Mirr = 0.85 and the remaining CMC parameters 
are chosen to be i?+ = 100, if = 0.1, and C= 1.0086. 

First, we compare initial data sets for an un- 
boosted black hole and for a boosted black hole with 
PJ7inaxA/iir = 1.77. Fig. [TO] shows the coordinate loca- 
tions of both the apparent horizon and the minimal sur- 
face for these two cases. The apparent horizon remains 
an approximate coordinate sphere, although its coordi- 
nate radius is reduced (recall that Mi„ is identical for 
the boosted and unboosted data set, which was achieved 
by reducing i?nis for the boosted case). Furthermore, the 
apparent horizon is offset from the excision sphere in a 
direction opposite to the boost P*, analogous to the be- 
havior of asymptotically flat inversion symmetric Bowen- 
York initial data fH. 
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FIG. 10: Coordinate locations of the apparent horizons and 
minimal surfaces, cut through the x-z plane, for a non- 
spinning, unboosted black hole (solid black circles), and for 
a non-spinning black hole boosted in the z-direction (dashed 
red circles). 
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FIG. 11: 2D Ricci scalar on the apparent horizon surface of 
a black hole with nmaxf'/Miir = 1.77. In this view, a wedge- 
shaped region has been removed from the front. 



To investigate the intrinsic geometry of the apparent 
horizon for the boosted black hole, we compute the Ricci 
scalar ^'^'^nM^^ fr om the 2-metric induced on the appar- 
ent horizon. Fig. [TT]plots this quantity; it is axisymmet- 
ric (as it must be), is maximum at the poles along the 
z-axis and minimum along the equatorial region. (Recall 
that (^^T^M^j. = 0.5 for a spherical geometry.) The nu- 
merically computed spin of this black hole is indeed zero, 
to machine precision. 

Fig. [12] shows the minimum and maximum of ^TZM?^^ 
as the boost parameter is varied in the range < 
f^max^'Z-^irr ^ 6.87. The minimum and maximum values 
of ^7?.inaxAf{J,^ when r^max-P/Mrr = 1-77 agrcc with those 
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FIG. 12: Dimensionless intrinsic geometry of the apparent 
liorizon (2D Ricci scalar, times the irreducible mass squared) 
versus Pfimax/Mirr for a single non-spinning black hole, with 
the ratio Rah/Ruis kept fairly constant. Shown are the max- 
imum and minimum immerical values computed on CMC 
slices. 



shown in Fig. [TT] 



D. Binary black hole initial data 

To demonstrate the generality of the approach that we 
have presented, we shall construct initial data for two 
black holes with mass-ratio approximately 2 : 1 and non- 
zero, arbitrarily oriented Bowcn-York spin and boost pa- 
rameters. First we describe how wc obtain input parame- 
ters for the elliptic solver corresponding to our particular 
physical parameters. First, we choose 7 = rms/rn = 0.8, 
which singles out a particular line of constant 7 in Fig. 
m for each black hole. Next, we pick KM for each 
black hole so that (i) its minimal surface is at least par- 
tially down the throat of the trumpet, which is near the 
turnover of the 7 = 0.8 curve and (ii) Eqs. (jlS]) and 
PS)) hold, i.e. before the turnover. With these criteria in 
mind, we choose K (a global parameter) to be 0.05 and 
the masses of black holes A and B to be, respectively. 
Ma = 2/3 and Mb = 1/3. Finally, from Eq. dM]) and 
Fig. El we find Rn,s/R+ = 8.1 x 10""* for hole A and 
Rn,s/R+ = 4.1 X 10-4 for hole B. 

We fix the overall length scale by setting i?+ = 300. 
This places the excision radii at Rms = 0.244 and 0.122 
for holes A and B, respectively, and the apparent horizon 
radu Rah ~ 1 (because from Fig. [H i?i„s/i?AH ~ 0.2). 
The coordinate locations of the two holes are then cho- 
sen to be {xa.Va.za) = (10,0,0) and {xb.Vb.zb) = 
(—20, 0, 0), and the center of mass of the holes is at the 
origin of the coordinate system. 

We take the spins to be Sa = {0,0, S a) and Sg ~ 
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FIG. 13: Convergence of the elliptic solver for the unequal 
mass binary black hole example shown in Fig. 1141 Shown is 
the volume L2-norm of the residual of f2 as A'" is increased, 
where A'^ is the cube root of the total number of collocation 
points. 



(53,0,0). Since we are adding significant spins, it is 
necessary to set CcS (defined in Eq. ([54]) ) equal to gcM^ 
for each hole, giving 0.569 for hole A and 0.142 for hole B. 
We take Sa = 0.4 and Sb = 0.1, then giving Ca ^ 0.0613 
and Cb = 0.0128. 

We take the boost parameters of the two black holes to 
be equal and opposite in the j/-direction, with magnitude 
Pa = Pb — 0.067. This gives approximate speeds of 
VA = ^^^^Pa/Ma = 0.24 and vb = n^^^PB/MB = 
0.48, with Ma = 2/3 and Mb = 1/3 as given above. 
With Ca and Cb not equal to zero, and not co-linear 
with 5", inversion symmetry is not possible (refer to the 
Appendix). Thus, we set = Qg = 0. 

Fig. [Tni shows exponential convergence of the volume 
L2-norm of the residual for the solution of the elliptic 
solver in this example, as the resolution of the numerical 
grid is increased. In addition, we have calculated the 
irreducible masses of the two holes and find values of 
0.53 for hole A and 0.27 for hole B. This gives a mass 
ratio of 1.96. 

Fig. [HI shows the conformal factor on the full compu- 
tational domain for the mass ratio 2 : 1 boosted, spinning 
binary black holes described above. The dark blue color 
at the outer edge shows that = at null infinity (to 
machine precision). In the middle, there is a prong- like 
feature, the tips of which arc the two black holes. It is 
evident that the conformal factor becomes quite small in 
the vicinity of the two black holes. 

Since our calculation of input parameters assumes 
spherical symmetry when our holes in fact have apprecia- 
ble spins and boosts, one expects the irreducible masses 
to differ somewhat from the values used for calculating 
the input parameters. This is indeed what we find (0.53 
vs. 0.67 for hole A and 0.27 vs. 0.33 for hole B). Finally, 
we find that the intrinsic geometry of each hole is dis- 
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FIG. 14: Conformal factor SI for a spinning, boosted binary 
black hole, with mass ratio 2:1. The inner boundaries are 
the minimal surfaces of the two black holes, and the outer 
boundary is null infinity. The maximum of Q is red, and 
equals 2.38. The minimum is dark blue, which is zero to 
machine precision. 



torted by the same amount. In particular, the minimum 
and maximum values of (^^7?.M^j. are, respectively, 0.37 
and 0.54 for each hole. 



IV. DISCUSSION 

In this paper, we have considered the conformal 
method on CMC hyperboloidal slices, focusing on gener- 
alizing the traditional Bowen-York data. There are two 
key aspects that make Bowen-York data easy to con- 
struct. First, for constant mean curvature (no matter 
whether A' = or i^T 7^ 0) , the momentum constraint de- 
couples from the Hamiltonian constraint, and the former 
can be solved first. Second, with conformal flatness, the 
momentum constraint simplifies to such an extent that 
analytical solutions are known: the symmetric, tracefree, 
divergence-free tensors with appropriate radial fall-off. 
Interestingly, the second aspect carries over from zero 
mean curvature to non-zero mean curvature, with the 
conformal factor now playing the dual role of turning the 
Hamiltonian constraint into an elliptic equation as well as 
compactifying The general Bowen-York conformal 

traceless extrinsic curvature still solves the momentum 
constraint analytically, where the only change necessary 
is addition of a spherically symmetric divergence- free ten- 
sor [the first term in Eq. ([M]) ]. 

Compared to the zero mean curvature case as usu- 
ally formulated, the hyperboloidal initial value problem 
has more free parameters, most notably the constants K 
(the mean curvature) and C (the coefficient of the spher- 
ically symmetric contribution to the conformal traceless 
extrinsic curvature), though a non-zero C is also consis- 
tent with zero mean curvature. Both these constants, as 
well as the minimal surface coordinate radius, have to 
be chosen carefully, and a significant portion of Sec. |TT] 
is devoted to working out permissible choices, and their 



consequences on the initial data under construction. 

As in the zero mean curvature case, hyperboloidal 
Bowen-York initial data trivially extends to multiple 
black holes with different spin- and boost-parameters for 
each black hole. Once again, one must be careful to 
choose the constants C (one for each black hole) and 
the radii of the excision boundaries, and Sec. Ill El gives 
simple rules how to do this. However, it is worth noting 
one significant difference. A single black hole must be 
centered at the origin of the conformal coordinates on a 
CMC hypersurface to be precisely Schwarzschild, since 
the outer boundary condition is imposed at a finite coor- 
dinate radius. A displaced black hole is not spherically 
symmetric. 

For hyperboloidal slices, the elliptic equation for the 
Hamiltonian constraint, Eq. pOp . is singular at the outer 
boundary J^^, where 17 0. Perhaps surprisingly, we 
have not encountered any difhculties when numerically 
solving this equation, for either single or binary black 
hole initial data. This is without any attempt to iso- 
late and explicitly cancel the singular terms in the equa- 
tion at future null infinity, as was advocated in [l8| . 
We suspect that the absence of numerical difficulties 
is related to the simple Dirichlct boundary condition 



0, and to the fact that the singular terms force 



the solution to also satisfy the von Neumann condition 
{dn/dR)y+ = -K/i, which follows from Eq. (HH), im- 
plying spherical symmetry to at least first order in an 
expansion away from null infinity. The freedom in the so- 
lution at the outer boundary necessary to accommodate 
a global solution of the elliptic equation also satisfying 
an inner boundary condition at _Rnis resides in a higher 
order term in the expansion of f2 away from Our 
spectral code never evaluates the Hamiltonian constraint 
Eq. Uni) right at J^+. Rinne has also had no dif- 
ficulty in solving the same elliptic equation with a finite 
difference code, as part of a constrained evolution scheme 
on CMC hypcrsurfaces [21 1. 

The Hamiltonian constraint equation is also singular 
at the inner boundary in the special case of a trumpet, 
for which = at i? = 0. This is a more challenging nu- 
merical problem, as discussed in Sec. Ill Fl Singular terms 
include some inside the Laplacian operator. Their can- 
cellation again uniquely determines the normal derivative 
of f2 there, but now that will have angular dependence 
if the spin is non-zero. The solution does not have a 
simple expansion in integral powers of R at the bound- 
ary, which makes it more of a challenge for our spectral 
methods to have the accuracy required to deal with the 
singularities in the equation. Still, we were able to ap- 
proach very close to the trumpet limit, at least in the 
spherically symmetric case, by using the analytic solu- 
tion for f2/i? at the boundary to formulate the boundary 
condition as a Dirichlet condition on f2 at a small, but 
non-zero, R and by adding extra collocation points near 
the boundary. If need be, reformulating the Hamiltonian 
constraint as an equation for fi/i?, as in Eq. (j47p . in a 
domain near the inner boundary and explicitly cancel- 
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ing the singular terms at the boundary should make it 
possible to manage exact trumpet boundary conditions 
numerically. 

One interesting aspect of hyperboloidal Bowen-York 
data lies in the physical interpretation of the spin pa- 
rameter 5' and the boost parameter P\ On asymptot- 
ically flat hypersurfaces (i.e. K = 0), one can evaluate 
the ADM integrals and find that the Bowen-York pa- 
rameter agrees with the ADM linear momentum, and 
that agrees with the ADM angular momentum. For 
hyperboloidal slices, the ADM formulas are not appli- 
cable. Nevertheless, a single unboosted spinning black 
hole, because it is axisymmetric, has a well-defined an- 
gular momentum which agrees with the spin parameter 
S"* (see Eq. (|43p ). The relationship of the boost parame- 
ters to the linear momentum is less clear. The conformal 
compactification leaves a rcscaling freedom R — > rjR un- 
specified (see Eq. ([24])). and as argued in Sec. Ill D"2| the 
boost parameter rescales as 77^ ^P', so that the vec- 

tor P* by itself has no physical meaning. However, one 
can define a scale invariant quantity, r^max^** , which may 
be considered a "quasi-local" linear momentum, at least 
when K is small. The proper interpretation of boosts on 
CMC hypersurfaces requires further analysis. 

When angular momentum is defined, we can consider 
the question of how large spins can be constructed with 
hyperboloidal Bowen-York data. In Sec. IIIIBl we have 
considered a sequence of black hole initial data with in- 
creasing spins, and Fig. [9] shows that black holes have 
been constructed with spin-extremality parameter'^ C = 
5/(2Mj2^J w 0.78. In contrast, standard Bowen-York 
data for a sin gle spinning black hole allows ^ < 0.83 
(Fig. 2 of Ref. [53 ), whereas conformally flat conformal 
thin sandwich data allows C, ^ 0.56 along the easily ac- 
cessible lower branch of solutions, and ^ ^ 0.87 along 
the upper branch (Fig. 7 of Ref. ^^)- We thus see that 
hyperboloidal Bowen-York initial data allows similarly 
large spins as the standard Bowen-York initial data (this 
includes the widely used puncture initial data psj as a 
special case). We have not tested the sensitivity of the 
maximum achievable ^ to variations of the other Bowen- 
York parameters C, K, of spinning black holes, but 
do not expect it to be large as long as KMi^i- is reasonably 
small. 

The simplifying assumptions of Bowen-York initial 
data appear to limit the ability to push towards near- 
extremal spins C w 1. To construct larger spins, one 
would have to give up these simplifying assumptions, 
most notably conformal flatness. An approach based on 
the extended conformal thin sandwich (XCTS) equations 
similar to Ref. seems very promising. Note that for 
the Schwarzschild spacetime, the space-time metric can 



^ We avoid the more widely used spin measure x = S/M'^, with M 
the Christoudoulou mass, because X ^ 1 due to the definition of 
the Christoudoulou mass, and because the Christoudoulou mass 
only has physical meaning for Kerr black holes. 



be conformally rescaled, resulting in conformal lapse and 
shift functions which are flnite at ^+ (sec Eq. ([2T|)). 
Thus, it seems quite likely that the XCTS equations 
rewritten in suitably rescaled variables can be used to 
construct more sophisticated hyperboloidal initial data. 
The XCTS-approach has another interesting feature. In 
this approach, the spins and boosts of the black holes are 
implemented by boundary conditions at the black hole 
horizons [1^, [lo, [Ml; the region of the initial data hy- 
persurface close to the black holes should only be mildly 
affected by the "warping up" of the CMC hypersurface 
at large radii as it approaches J^~^ . Therefore, within the 
XCTS framework, it might be easier to interpret a boost. 
This will be a topic of future research. 
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APPENDIX: THE INVERSION SYMMETRY OF 
THE CMC INITIAL VALUE PROBLEM 

In discussions of the conformally flat initial value prob- 
lem, there has been considerable interest in inversion 
symmetric initial value data (see (6^ and references 
therein). The issue of inversion symmetry arises when 
the initial hypersurface contains a minimal surface at a 
conformal radius R = Rms- As R decreases below this 
value, the physical radius r increases and becomes in- 
finite in the limit R —^ 0. The hypersurface may or 
may not be symmetric under the inversion transforma- 
tion R — !• R^g/ R. Imposing inversion symmetry on the 
initial data, and requiring that it be preserved during 
subsequent evolution, can lead to relatively simple exci- 
sion boundary conditions at P = . Inversion symme- 
try was discussed in the original Bowen and York paper 
[22{ . and has been exploited in much of the numerical 
work based on the Bowen-York class of solutions to the 
initial value problem on maximal hypersurfaces. In this 
Appendix, we show that the conditions on the solution for 
the conformal tracelcss extrinsic curvature tensor which 
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lead to inversion symmetry are the same on CMC hyper- 
surfaces as they are on maximal hypersurfaces, noting, 
however, that the most general Bowen-York solution of 
the conformal momentum constraint equation does not 
admit inversion symmetry. 

We start with the Hamiltonian constraint equation in 
the form given in Eq. (|47p as an equation for the scale- 
invariant variable U = fl/R. A rearrangement of terms 
gives a form in which the possibility of inversion symme- 
try is manifest: 



3 

2U 



dU" 



1 



(A.l) 



Under an inversion transformation, the first two terms 
on the right-hand side of Eq. (|A.2p transform into each 
other provided that = ±i?^jgP\ The next three 
terms do not depend on R and are therefore trivially 
inversion-symmetric. Symmetry of the second to last 
square-bracket requires Q' = +i?j^gP% while symmetry 
of the last square-bracket requires = —R^^P'^ unless it 
vanishes because the boost and spin vectors are co- linear. 
Without any restrictions on the Bowen-York parameters 
(C, P*, 5') , there is no choice of the which guarantees 
an inversion-symmetric W and therefore no guarantee of 
an inversion-symmetric solution of the Hamiltonian con- 
straint equation. This result does not depend on the 
value of if . If we set C = 0, we recover the inversion 
symmetry result as usually stated for maximal hyper- 
surfaces, that the "minus" form of inversion symmetry 
applies for general spin and boost. 



where A is the Laplacian operator and V is the gradient 
operator on the unit two-sphere. The only term not ob- 
viously symmetric under the inversion transformation is 
the term involving the source term W. If and only if W 
is inversion symmetric, W{R,9,ip) = W (R'^^^/R,0,(j)), 
will the solution for U, subject to the minimal surface 
condition dU/dR = at i? = Pms, be inversion symmet- 
Vic, U{R,0,^)^U{RlJR,0,cP). 

The generalized Bowen-York solution for Aij is given 
in Eq. ^9^. From this, we find 



-p2 [P^Pk + 2 (P^O (P'n,)] 



-P-2 ^QkQ^ + 2 (g^n,) (Q%,)] 

- 9 [P'^Qk - 4 (P'n,) (Q%,; 
18 (e,,fc5%^-) (e"""^™n„) 
1 
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RiP^ni,] 



R 
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R{e,,,P'S^n'') 



R 



(A.2) 
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